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ABSTRACT 

We use Smoothed Particle Hydrodynamic simulations of cold, uniform density, 
self-gravitating filaments, to investigate their longitudinal collapse timescales; these 
timescales are important because they determine the time available for a filament to 
fragment into cores. A filament is initially characterised by its line-mass, /i^, its radius, 
Rq (or equivalently its density = pt^/TiR^^)^ and its aspect ratio, (= Z^jR^^ 
where is its half-length). The gas is only allowed to contract longitud inally, i.e. 
parallel to the symmetry axis of the filament (the 2 ;-axis). fPon et all (j2Q12f ) have con¬ 
sidered the global dynamics of such filaments analytically. They conclude that short 
filaments {A^ < 5) collapse along the 2 ;-axis more-or-less homologously, on a time- 
scale ^ 0.44 Aq (GPq)“^/^; in contrast, longer filaments {A^ >5) undergo end- 
dominated collapse, i.e. two dense clumps form at the ends of the filament and converge 
on the centre sweeping up mass as they go, on a time-scale ^ 0.98 (Gp^ 

Our simulations do not corroborate these predictions. First, for all A^ >2, the collapse 
time satisfies a single equation 

{0A9 + 0.26A^){Gp^)-y\ 

which for large A^ is much longer than the Pon et al. prediction. Second, for all A^ >2, 
the collapse is end-dominated. Third, before being swept up, the gas immediately 
ahead of an end-clump is actually accelerated outwards by the gravitational attraction 
of the approaching clump, resulting in a significant ram pressure. For high aspect ratio 
filaments the end-clumps approach an asymptotic inward speed, due to the fact that 
they are doing work both accelerating and compressing the gas they sweep up. Pon et 
al. appear to have neglected the outward acceleration and its consequences. 

Key words: ISM: clouds - ISM: kinematics and dynamics - ISM: structure - stars: 
formation 


1 INTRODUCTION 


The u biquity of non-spherical s t ructu re s in molecular clouds 
(e.g. ISchneider Elmegreen 1 19791 : [Lada. Alves Lad^ 


1 19991 : iMverd ~ 20091 : lAndre et al.l l201Ql i lead to 

many papers on t he effe cts _ of dimensionality on 

global collapse (e.g. Bastien 19831: iBastien et al. | l991 : 

hiboii 


Burker^jMlartmann|200jjPQnJj]^n^tone jMlei^^ 


Toala. Vazquez-Semadeni G6mejl2012l : |Pon et aDl2012l k 


It has been shown analytically that non-spherical struc¬ 
tures collapse on lon ger timescales than spherical clouds of 
the s ame density (e.g. Burkert Hartmann|[20o3 : |Pon et al.l 
2012lk and this has been corroborated by si mulations (e.g. 
BastienI [l 9831 : Ivizauez- Semadeni et aLllioO^ k In cylindrical 
clouds - hereafter filaments - there are two possible longi- 
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tudinal collapse modes (i.e. collapse modes parallel to the 
symmetry axis of the filament): a filament may collapse ap¬ 
proximately homologously (hereafter homologous collapse), 
or clumps may form at the ends and move inwards sweeping 
up material as they go (hereafter end-dominated collapse; 

_ 

Recently, I Pon et aT ] il2ni2h have explored analytically 
the effect that these two modes have on the freefall collapse 
time for a uniform filament. During homologous collapse the 
density of a filament remains approximately uniform but 
changes with time, so all points within t he filament have 
approximately the same collapse timescale. IPon et al.l (l2012l) 
calculate that the homologous collapse timescale is 

^HOM - O.UA^ • (1) 

Where is the initial volume density of the filament, Aq 
is the initial aspect ratio defined as Aq = Zq / R ^, Zq and 
Rq are the filament’s half-length and radius respectively. 
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Figure 1. (a) The volume-density profile, and (b) the velocity profile, along the 2 ;-axis at time t = 0.55 Myr, from the simulations of 
filaments having = 10 Mq pc“^, = 0.1 pc and = 3 (dashed curves) or = 10 (solid curves). The signatures of end-dominated 

collapse are clearly seen in both, viz. strong density peaks and supersonic inwards motions at the extremes, representing the end-clumps 
that dominate the dynamics. The horizontal dotted line in figure (a) is the initial volume density; interior material is relatively unaffected 
by the collapse until the density peaks reach it. 


The exact equation describing the longitudinal accelera¬ 
tion inside a filament is non-linear, and the non-linear terms 
become very important at the ends of a filament. Conse¬ 
quently the ends of the filament are accelerated to much 
greater speeds than in homologous collapse, cau sing them 
to swe ep up interior matter as they overtake it. |Pon et ^ 
(|2012l ) calculate that in this case the collapse time is 

c. 0.98 4/qGpo)-'/^ (2) 

Since the collapse times for t he two modes depend on 
Aq in different ways, |Pon et al ] (I 2012 I) conclude that fil¬ 
aments collapse on whichever is the shorter timescale, i.e. 
short filaments collapse approximately homologously, on a 
time-scale given by Eqn. o, and long ones collapse via 
the end-dominated mode, on a timescale given by Eqn. 
Comparing Eqns. o and they infer that the switch 
occurs at Aq ~ 5. 

The two modes should be distinguishable by their kine¬ 
matics. A short filament collapsing homologously has an ap¬ 
proximately linear velocity profile along its long axis, while 
a long filament collapsing via the end-mode has low inward 
velocities except at its ends, where the inward velocity in¬ 
creases dramatically. 

In this paper we present numerical simulations of ini¬ 
tially uniform filaments to evaluate their freefall time-scales, 
and to inv e stigat e the collapse-mode dichotomy presented in 
|Pon et all (|2012h . In 32] we detail the simulation setup and 
the initial conditions used; in 32] we present the results of 
these simulations; in SI we discuss their significance and 
support the discussion with a semi-analytical model of the 
system dynamics; and in 33 we summarise our conclusions. 


2 NUMERICAL SETUP 

The simulations presented in this paper are performed using 
the Smoothed Particle Hydrodynamics (SPH) code GAN- 
DALE (Rubber et al. in prep). The simulations invoke both 


gravity and hydrodyna mics, with an isotherma l equation 
of state. Grad-h SPH (|Price Monagfi^ l2004h is imple¬ 
mented, with ?7 = 1.2 0 , so that a typical particle has ~ 56 
neighbours. Although we are interested in freefall collapse, 
the gas is given a non-zero temperature, T = IK, in order 
to capture shocks, and to stop artificial clumping of SPH 
particles. Eor molecular gas with solar composition, this re¬ 
sults in an isothermal sound speed of 0.07 kms“^. Running 
identical simulations at different low temperatures produces 
only small variations in the collapse times, so we are satisfied 
that T = 1 K gives a reasonable approximation to free-fall. 

The code is modified so that radial motions are sup¬ 
pressed, i.e. Vx = Vy = 0 for all particles. Thus, the longitu¬ 
dinal collapse of a cylindrical cloud with super-critical line 
mass can be simulated without it collapsing radially. Such 
a constraint has no effect on the longitudinal motion of the 
gas and allows the simulati ons to conform to the constant 
radius condition adopted by |Pon et all (|2012l ) . 

A hexagonal close-packed grid of particles is used to 
produce a uniform density. The grid is rotated around the 3 
axes by random angles to ensure that the grid’s symmetry 
axes are not aligned with the Gartesian axes, and a cylindri¬ 
cal filament is then cut from this grid. The initial conditions 
are characterised by 

Aq, the aspect ratio of the filament; 

/Xq, the mass per unit length of the filament; and 

Rq, the radius of the filament. 

Aq is varied between 1 and 20. A filament is defined 
as a structure with a n aspect ratio of at least 3 (e.g. 
IPanopoulou et aDl2014l ): aspect ratios below this value are 
only simulated for completeness. /Xq is varied between 
2 MqPc“^ and bOM^pc”^. is varied between 0.05pc 
and 0.15 pc. This results in a range of initial number densi- 


^ ?7 is a constant used to calcula te smoothing lengt hs self- 
consistently in Grad-h SPH (see e.g. iHnbber et al .11201 lib It de¬ 
termines the average number of neighbours. 
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Figure 2. The filled circles show values of the collapse time, t^oL ’ 
obtained when (a) the aspect ratio, , is varied and and Rq 
are fixed at their fiducial values; (b) the line-density, , is varied 
and Aq and are fixed at their high aspect ratio fiducial values; 
(c) the radius, R ^, is varied and and /j,^ are fixed at their 
high aspect ratio fiducial values. The high aspect ratio fiducial 
values are A^ =10, =10Mq pc“^, and R^ =0.1 pc. The blue 

continuous lines are the predictions of Eqn. ([3] ). The red dotted 
and green dashed lines are the predictions from IPon et al ] il2ni2l~i 
(Eqns. [T] and [2] respectively) 



Figure 3. The filled circles show values of the collapse time, ^col ? 
obtained with different aspect ratios, Aq , when the radius, Rq , is 
fixed at 0.1 pc, and the line-density is fixed at /Xq = 35Mq pc“^ 
or [Iq = 50 Mq pc“^. The continuous lines are the predictions of 
Eqn. [3] with = p^/ttR’^. 

ties, from 5 x 10^ cm“^ to 1 x 10^ cm“^, for molecular gas 
of solar composition. 

Much of our discussion will consider two fiducial cases: 
a high aspect ratio case characterised by A^ = 10, /Xq = 
10 Mq pc“^ and R^ = 0.1 pc and a low aspect ratio case 
characterised by A^ =3, /Xq = 10Mq pc“^ and R^ =0.1 pc. 

The simulations are all run with the same initial 
volume-density of SPH particles, ~ 800,000 pc“^, 

to ensure that they have comparable spatial resolution, 
~ 0.01 pc; the mass resolution therefore varies between 
0.008 Mq and 0.20 M^. Having comparable spatial resolu¬ 
tion (rather than comparable mass resolution) is appropriate 
since we are concerned here with the longitudinal density- 
and velocity-profiles (rather than fragmentation into con¬ 
densations). To check convergence, we have repeated three 
simulations with rXgp^ ~ 2,400, 000 pc“^; the resulting col¬ 
lapse times, density- and velocity-profiles differ by < 1%. 
We therefore believe that rigp^ ^800, 000 pc“^ delivers suf¬ 
ficiently high resolution. We have also simulated the high 
aspect fiducial case (A^ = 10, /Xq = 10 pc~^, R^^ — Q -1 Pc) 

using a different SPH code fSERENi lHubber et al.ll201lh . and 
obtain almost identical results, so we are confident in the fi¬ 
delity of the simulations. 


3 RESULTS 

For all aspect ratios, 2 < Aq < 20, all line-densities, 
2 MqPc“^ ^ Mo ^ 50MqPc“^, and all radii, 0.05pc < 
Rq ^ 0.15 pc, the global collapse of the filament is end- 
dominated. The ends experience the highest acceleration, 
and this leads to the formation of end-clumps which then 
converge on the centre, sweeping up material as they go. 
Fig. □ shows the density and velocity profiles along the z- 
axis, for both the high and low aspect ratio fiducial cases. 
The dense end-clumps and their supersonic inward velocities 
are clearly seen. 

In all the simulations, we estimate the collapse time, 
tcoLJ t)y identifying the moment when the two converging 
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Figure 4. (a) The velocity of the z > 0 end-clump as a function of time for filaments with initial line-density = 10 Mq pc“^, radius 
Rq = 0.1 pc, and different initial aspect ratios, = 3, 5, 10, and 20. In the high aspect ratio simulations the end-clump accelerates 
inwards for the first ~ 1 Myr, but then approaches a terminal velocity, which is only weakly dependent on . The sudden increase in 
the negative velocity at the end of each profile is due to the two end-clumps becoming so close that the motion is dominated by their 
mutual gravitational attraction; there is essentially no gas left in between them to sweep up. (b) The total kinetic energy and the work 
done on the gas plotted against time for both the high and the low aspect ratio fiducial cases (p,Q = 10 Mqpc“^, Rq = 0.1 pc and 
Aq = 10 or Aq = 3). For the Aq = 10 case, the end-clumps approach their terminal velocity, around 1 Myr; thereafter, the rate at 
which work is done exceeds the rate of increase of the kinetic energy. This is not seen in the Aq = 3 case due to the filament collapsing 
before terminal velocity is reached; the abrupt shift in energies at ~ 1.1 Myr is due to the end-clumps meeting at 2 ; = 0. 


end-clumps merge, i.e. when there is a single central density 
maximum at z = 0. The filled circles on Fig. [2^ show values 
of tcoL obtained with different Aq, when /Xq and R^ are 
fixed at their fiducial values. The filled circles on Fig. [JJo 
show values of obtained with different /Xq, when Aq 

and Rq are fixed at their high aspect ratio fiducial values. 
The filled circles on Fig. ^ show values of ^qql obtained 
with different R^, when Aq and /Xq are fixed at their high 
aspect ratio fiducial values. The blue continuous curves show 
that all these results are well fitted with 

icoL = (0.49 + 0.26^o) . (3) 

where = /Xq /ttR^ — and hence they are n ot well fitted 
with the relations derived by |Pon et all (|2012l l (i.e. Eqn. [1] 
for Aq < 5, and Eqn. for Aq > 5). 

Eigure [3] compares the collapse times obtained from 
simulations with R^ = 0.1 pc, /Xq = 35 Mq pc“^ or /Xq = 
50MqPc“^, and a range of values for Aq (filled circles) 
with the predictions of Eqn. 0 continuous green and blue 
lines). Again the agreement is very good. 


4 DISCUSSION 

The gravitational acceleration felt by particles off the sym¬ 
metry axis is lower than that felt by those on the axis, an 
effect we discuss below in an This results in a curved 
end-clump, the portion of the end-clump on the axis being 
roughly 0.01 pc ahead of the end-clump off axis. To deter¬ 
mine the velocity of the end-clump, we locate the highest- 
density SPH particle that is at least 0.04 pc off axis, and 
then average the velocity of all the SPH particles within 
z zb 0.02 pc of this particle. By selecting a SPH particle off 
axis we are able to avoid the atypical central region and lo¬ 
cate the end-clump more accurately. The resulting velocity 


is insensitive to the precise choice of zb 0.02 pc and of 0.04 pc 
as the off axis requirement. 

Eig. [3^ shows the velocity of the 2 : > 0 end-clump as 
a function of time, for the four filaments with initial line- 
density /Xq = 10 MqPc“^, radius R^ = 0.1 pc and initial 
aspect ratios, Aq = 3, 5, 10, and 20. The velocity profiles 
are roughly independent of aspect ratio. Eor large values of 
Aq a terminal velocity is reached before the filament has 
time to fully collapse. 

A terminal velocity arises because, once an end-clump 
becomes sufficiently massive, it starts to dominate the grav¬ 
itational field experienced by the gas just ahead of it, and so 
this gas flows outwards to meet the end-clump (see Eig. [5^). 
This means that more work has to be done reversing the 
velocity of the swept-up gas, and the swept-up gas is more 
strongly compressed, so more work is done on it. The up¬ 
shot is that less gravitational potential energy gets converted 
into kinetic energy, basically just enough to accelerate the 
newly swept-up gas up to the terminal velocity. This is illus¬ 
trated on Eig.[4lD, where we plot, for the simulations with the 
high and low aspect ratio fiducial parameters ( R^ = 0.1 pc, 
/Xq = lOM 0 pc“^ and Aq = 3 or Aq = 10), the evolution 
of the net kinetic energy and the net work done. At early 
times, the kinetic energy grows fastest, but once the termi¬ 
nal velocity has been reached for the Aq = 10 case, around 
~ 1 Myr, work done takes over, and the kinetic energy in¬ 
creases roughly linearly with time, due to the roughly linear 
increase in the mass of the end-clumps as they sweep in at 
constant velocity. This large growth in work done is absent 
for the Aq = 3 case; this is to be expected because the 
end-clumps merge in the centre before terminal velocity is 
reached. However, the same dynamics are present, the end- 
clump gains enough mass to dominate the gravitational field 
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Figure 5. (a) The gas velocity along the 2 ;-axis at time t = 
1.2Myr, for the filament with Aq = 10, /Tq = 10Mqpc“^ and 
Rq = 0.1 pc. The gas just interior to the converging end-clumps 
is moving outwards, thereby increasing the ram-pressure that this 
gas delivers to the end-clump when it is swept up. (b) The density 
profile of the same filament at the same point in time. The density 
peaks (i.e the end-clumps) are at the same position as the velocity 
peaks. Note the uniformity of the interior density, it has remained 
roughly constant and uniform since the start of the simulation. 


in its vicinity and the gas just ahead of it begins to move 
outwards to meet it. 


4.1 Semi-Analytical Model 


As well as considering the system in terms of energy, we 
can also develop a semi-analytical model on the basis of 
the forces involved. A terminal velocity is reached, for high 
aspect ratio filaments, because the gravitational force on 
the end-clump is approximately balanced by the force of 
ram-pressure due to the end-clump sweeping up the outward 
moving interior material. The net acceleration of the end- 
clump is therefore 


(Pz 


- a{2TiG PqRq) + 


Po 

M(i) 


(4) 


The first term on the righthand side of Eqn. 0 is the 
gravitational acceleration that the end-clump experiences 
due to the interior material. {^ttGpqRq) is the acceleration 


Figure 6. The 2 ;-velocity of the gas just interior to the 2 : > 0 
end-clump as a function of time. At early times, this gas flows 
inward with the end-clump, but at later times, it is attracted 
gravitationally to the end-clump and moves outwards to meet it. 
Close to the end of the collapse, the interior material begins to 
be gravitationally attracted by the 2 ; < 0 end-clump resulting in 
the velocity maximum at ~ 2.2 Myr. 


derived by |Pon et aP (|2012l ) from consideration of a point 
on the symmetry axis of the filament, a is a correction factor 
that takes account of the fact that the end-clump extends 
from the symmetry axis, r = 0, to r^R^, where the gravi¬ 
tational acceleration is lower. 

The second term on the righthand side of Eqn. © is the 
acceleration experienced by the end-clump due to ram pres¬ 
sure; is the relative velocity between the end-clump 

and the material about to be swept up, and M(t) is the mass 
of the end-clump at time t. If the material interior to the end- 
clump were stationary, we could put ^dzjdt, where 

dz/dt is the velocity of the end-clump at time t. However, 
the interior material is not stationary, as seen in figure [5^, 
and so we must put —dzidt — dz^jdt^ where dz^jdt is 

the velocity of the material just interior to the end-clump at 
time t. Equ ation m is equivale nt to the momentum equation 
used in the IPon et alJ (l2012h paper. It is their assumption 
that Ar'pp^L — dzidt (i.e the interior material is stationary) 
which distinguishes their analysis from ours. 

Eig. [6] shows that the velocity of the interior material 
is relatively low, and it moves inwards and then outwards. 
If one neglects the net motion of the interior material, then 
the mass of the end-clump is M{t) — z{t))^ and the 

acceleration of the end-clump becomes 




— a {2 7 t GPq Rq) + 


(dz/dt — dzj/dty 
(Zo - z(t)) 


(5) 


The assumption that the end clump’s mass is approx¬ 
imately the same whether the interior material moves or 
not is shown to be valid in Eig. [71 Here, the mass of the 
z > 0 simulated end clump is compared with the function 
M(t) = Pq{Zq — z{t)) for the high aspect fiducial case. The 
masses agree with each other within ~ 10%. This is due to 
the interior density staying roughly constant during the col¬ 
lapse. At such low temperatures the end-clumps are moving 
supersonically and thus no density waves disturb the mate¬ 
rial interior to the clumps. Moreover, most of the interior 
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material is unaffected by the gravitational pull of the end 
clumps and is approximately stationary until an end-clump 
approaches (Fig. [5]). 

To solve Eqn. for z{t) (the location of the end- 
clump), one must first find dz-^jdt, the speed of the material 
about to be swept up by the end-clump. The acceleration at 
Zj is 

^ - -a(27rGpo) |2«j(i) + (fig + 

- {rI + {z{t)+zdt)fy^^] 
+ /3(27rGpo)(^o-2W)|(^2 + 

{z{t)-zdt)) ] 

{Rl + {z{t)-zW)ry^^ i 

(6) 

The first term on the right hand side of Eqn. dni the 
one preceded by a) is the contribution to the gravitational 
acceleration from the material between the end-clumps, 
i.e. the acceleration at a point z^(t) in a uniform-density 
cylinder of radius betw een ± 2 :(t), where z^{t) < z{t) 
(|Burkert Hartmann l20Q4l i . We have again neglected mo¬ 
tion of this material, and included the factor a to take ac¬ 
count of the fact that the material away from the symmetry 
axis experiences a lower acceleration than the material on 
the axis. 

The second term on the right hand side of Eqn. dni the 
one preceded by (3) is the contribution to the gravitational 
acceleration from the end-clumps. This term is obtained by 
treating the end-clumps as discs of uniform surface-density, 
considering a point on the symmetry axis at z^ {t ), and then 
including a factor /3 to allow for the fact that material away 
from the symmetry axis experiences a lower acceleration 
than the material on the axis. 

Eqns. © and are second order non-linear cou¬ 
pled differential equations. We have solved them numer¬ 
ically for the high aspect ratio fiducial case — 10, 

/Xq = 10Mq pc“^, Rq =0.1 pc), using the first order Euler- 
Cromer method and 5t = 0.003, 0.010, and 0.020 Myr; the 
results obtained with these different 5t values are essentially 
indistinguishable, indicating that the Euler-Cromer method 
is converged. The interior of the filament is divided into 500 
evenly spaced Lagrangian discs and the velocity of each disc 
calculated at each time-step by integrating Eqn. Q. The 
velocity of the disc just interior to the end-clump is then 
used to solve Eqn. ©• The correction factor a can be eval¬ 
uated from the simulation data by averaging over the SPH 
particles in the end-clump, to obtain a = 0.60. The correc¬ 
tion factor {3 is obtained by averaging over the SPH particles 
making up discs of gas interior to the end-clumps, the first 
term of equation [6] is then subtracted. This results in a value 
of/3 = 0.4. 

Eig. [8] compares the results of the semi-analytical model 
with the simulation, for the high aspect ratio fiducial case 
{Aq — 10, /Xq = 10 Mq pc“^, Rq =0.1 pc) and the low aspect 
ratio fiducial case {A^ =3, /Xq = 10 Mq pc“^, R^ =0.1 pc); 
we have also made this comparison for the cases with Aq = 
10, Rq = O.lpc, and /Xq = 35MqPc“^ and SOM^pc”^, 
and the agreement is equally good. This suggests that the 
semi-analytical model is capturing the important physical 



Figure 7. The mass of the z > 0 simulated end-clump is com¬ 
pared with the end-clump mass expected if the interior material 
was stationary, M(t) = iIq{Zq — The masses are within 

~ 10% over the course of the simulation. 


processes, and that the factors a = 0.60 and /3 = 0.4 are valid 
correction terms. 

If Eqn. (|5)) accurately describes the acceleration of an 
end-clump, then, once the end-clump approaches its termi¬ 
nal velocity and its acceleration falls towards zero, we must 
have 


(dzjdt — dz-^jdt)^ 
(Zo - z(t)) 


a2 7T G Pq Rq . 


(7) 


Since the righthand side of this equation is constant, a plot 
of \og^^{dz/dt — dz-^/dt) against \og_^^{ZQ — z(t)) should 
have slope 0.5. This is shown in figure [9] for the high aspect 
ratio fiducial case {A^ = 10, /Xq = 10 pc“^, R^ =0.1 pc) 

between 1.2 Myr and 2.2 Myr. ft is well fit by a straight line 
with a gradient of 0.64, close to the expected 0.5. The de¬ 
viation is due to the fact that the end-clump is not quite 
moving at a constant velocity; it experiences a small accel¬ 
eration. There is a small oscillation about the straight line, 
indicating that, once the terminal velocity is approached, 
the dynamics is stable: if the end-clump is moving slightly 
slower than the terminal velocity, the acceleration due to 
ram pressure is slightly less than the gravitational accelera¬ 
tion, and it speeds up; conversely, if the end-clump is mov¬ 
ing slightly faster than the terminal velocity, the acceleration 
due to ram-pressure is slightly greater than the gravitational 
acceleration, and it slows down. 

Eqn. (|3]) agrees with the generic theorem that the free- 
fall collapse time for a uniform-density filament is longer 
than for a sphere of the same un iform dens i ty. Ho wever, it 
disagrees wi th the prediction s of jPon et ^ (120121 1 (Figure 
UK). Overall, |Pon et aT ](|2^ predicts somewhat longer col¬ 
lapse times for short filaments, A^ <9, an d shorter collapse 
times for longer filaments, A^ > 9. The |Pon et all (|2012l i 
collapse times for short filaments are too long, because their 
analysis does not take account o f the fact that ev en in this 
case there are end-clumps. The |Pon et ahl (|2012l ) collapse 
times for long filaments are too short, because their analysis 
does not take account of the interaction between an end- 
clump and the material it is about to sweep up, as analysed 
in SI We do not consider the discrepancies between our 
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Figure 8. A comparison between simulation data and our semi-analytic model for a filament characterised by = lOM^pc”^, 
Rq =0.1 pc and either = 10 or =3. (a) The position of the end-clump as a function of time. It starts at t = 0.2 Myr when the 
end-clumps are clearly established, (b) The velocity of the end-clump as a function of time. The velocity of A^ = 3 end clump has been 
shifted downwards by 0.1 kms“^ for clarity. Simulation and model profiles fit reasonably well for both initial aspect ratios. 


results and those found in |Pon et all (|2012h for the range 
Aq < 3 as such a cloud is outside our definition of a fila¬ 
ment. 


5 CONCLUSIONS 

|Pon et ^ ] (l2012ll have suggested that there are two modes 
for the freefall collapse of a uniform-density filament: short 
filaments, A^ < 5, collapse homologously, on a timescale 
that varies as A^ (see Eqn. [1]); longer filaments undergo 
end-dominated collapse, on a timescale that varies as 
(see Eqn.[2|. We find no such dichotomy. Rather, end effects 
are evident for all aspect ratios, and just one equation, 

^coL = (0.49 + 0.26 Aq) (Gpq) (|3)) 

predicts the free-fall collapse times obtained in our simula¬ 
tions to within ±3% for 2 < Aq < 20, 2Mq pc“^ ^ Mo ^ 
50 Mq pc“^, and 0.05pc < < 0.15pc. 

The linear relationship between large 

aspect ratios, is due to the fact that the end-clumps quickly 
approach a terminal velocity before collapse finishes. The 
same linear relationship still holds for small aspect ratio fil¬ 
aments despite the lack of a terminal velocity. Hence, the lin¬ 
ear relationship is ultimately a result of the dynamics of the 
system which are independent of initial aspect ratio. Before 
being swept up, the gas immediately ahead of an end-clump 
is accelerated outwards by the gravitational attraction of 
the approaching end-clump. This outward velocity increases 
the amount of ram pressure exerted on the end-clump, and 
the amount of work need e d to a ccelerate and compress the 
swept up gas. |Pon et aT ] (I2012I) appear to have neglected 
this outward moving gas and its consequences. 

We have constructed a semi-analytical model which de¬ 
scribes the acceleration of the gas in a collapsing filament, 
using only gravity and ram-pressure. This model is able to 
reproduce the end-clump position and velocity seen in the 
simulations and accurately describes the evolution of the 
end-clumps over a range of initial conditions. 



Figure 9. A test of the validity of equation [7| \og^^(dz/dt — 
dzjjdt) is plotted against \og^^{ZQ — z{t)) for the high aspect 
ratio fiducial case. The simulation data are well fit by a straight 
line with a gradient of 0.64, close to the expected 0.5. The oscil¬ 
lations around the linear fit show that once the terminal velocity 
is approached, the dynamics is stable. 
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